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Abstract. 

The connection between collimation and acceleration of magnetized rela- 
tivistic jets is discussed. The focus is on recent numerical simulations which 
shed light on some longstanding problems. 



1 Introduction 

Understanding the physics of relativistic jets, notably their acceleration and 
collimation, has been a major area of study for several decades. It is widely 
agreed that jets are a magnetic phenomenon. The basic paradigm involves a 
bundle of field lines that is attached at its base to a central compact star - black 
hole (BH) or neutron star (NS) - or to an accretion disk. Rotation of the star 
and/or disk causes a helical outgoing magnetohydro dynamic (MHD) wave which 
accelerates any gas that is frozen into the field lines. 

A gre at deal of analytical work has been done on jets , and the reader is 
refe rred tolBegelman k O dl994h.lyiahakis k Konkl d2003allbh. Beskin & Nokh- 
rina (|2006l ) . iNaravan et al.l (|2007l ) . iLvubarskvl (|2009al lbl) . and references therein, 
for an introduction to the vast literature. However, the equations describing 
MHD jets are nonlinear and the problem is complicated, so there is a limit to 
what can be accomplished purely analytically. Fortunately, in recent years, nu- 
merical investigations have begun to contribute to the field, facili t ated by the 



development of robust relativistic MHD codes (IKoide et al.l 120001: iKo missarov 



2001 



McK innev fc Gammiel 12004]; |Pe Villiers et al.ll2005l; iMcKinnev fc Naravan 



2007a Ibl: IKomissarov et al]l2007U2009l: iTchekhovskov et al.ll2008l. l2009aHbHcf) 



2 Force-Free Jets 

Before considering the full MHD problem, it is useful to focus first on force-free 
jets. The force-free approximation is a simplification of ideal MHD in which we 
include electric and magnetic fields, and the corresponding charges and currents, 
but we ignore the inertia of the gas. The problem thus reduces to pure electro- 
dynamics in a perfectly conducting medium. This very simple approximation 
provides surprisingly useful insight into the general problem. 
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2.1 Collimation 



Michel (jl973l ) derived an analytical solution for a force-free wind from a rotating 
star. He assumed a split monopole configuration in which the magnetic field is 
radial and has a constant strength at the surface of the star, with lines pointed 
into the star over one hemisphere and out over the other hemisphere. 

The rotation of the star causes the magnetic field to develop a toroidal 
component which dominates with increasing distance: = — {VtR/c)B p , 
where is the angular velocity of the star, R is cylindrical radius with respect to 
the rotation axis, and B p is the field strength in the poloidal (r0) plane. However, 
even though B^ dominates over B p at large radii {—B^jB p = ViR/c ^> 1), 
nevertheless there is no change in the poloidal geometry - the poloidal field 
continues to remain radial. This is surprising since the tension associated with 
toroidal field curvature, the so-called hoop stress, is very large. Why does this 
not collimate the jet? It turns out that hoop stress is pefectly canceled by an 
electric force in the opposite direction, and hence there is no net tendency to 
collimate the flow. (This is a purely relativistic phenomenon - the electric field 
is negligible in nonrelativistic MHD.) 

A sim i lar re sult i s seen in the parabo l oidal force- free solution derived by 
Blandfordl (|l976l ) and iBlandford Znaiekl (Il977f). and a more general class of 



self-similar force- free solutions studied bv lNaravan et al.l (|2007l ). Although hoop 



stress is very large at large distances from the base of the jet, once again it has 
small or no effect on the collimation of the jet. 

All this suggests that a relativistic jet is unlikely to self-collimate. It needs 
to be collimated by an external agency. Assuming that the jets in active galac- 
tic nuclei (AGN) and X-ray binaries (XRBs) are associated with field lines from 
a spinning accreting BH, these jets must be confined by a wind from the sur- 
rounding accretion disk. In the collapsar model of gamma-ray bursts (GRBs), 
the confinement is presumably due to the stellar envelope of the collapsing star. 
If no confining medium is present, as in the case of relativistic outflows from 
radio pulsars, the flow is nearly radial. 

2.2 Acceleration 

Force- free jets accelerate very efficiently In the simplest geometries, viz., split 
monopole and paraboloidal, the Lorentz factor of the outflow varies asymptot- 
ically as 7 ~ £lR/c, increasing linearly with cylindrical radius. In the case of a 
paraboloidal jet, 7 oc z 1 ! 2 (since R oc z 1 / 2 ), where z is distance along the jet. 
For more general geometries such as R oc z^ 2 ~ v ^ 2 ( Naravan et al. 20071 ). 7 



depends both on the value of QR/c and on the poloidal curvature of field lines. 
As a result, the run of 7 with z exhibits two distinct asymptotic regimes (see 
iTchekhovskov et al. 2008 for details). Nevertheless, at any given z, the field line 
with the largest Lorentz factor still satisfies 7 oc z 1 / 2 . Thus, this simple and 
convenient scaling appears to be fairly general. 

2.3 Jet Power 

Since a force-free jet is purely electromagnetic, the energy flux is given by the 
Poynting vector E x B/4ir (we set c = 1). For a rotating axisymmetric jet, 
E = (VLR/c)B p and is pointed perpendicular to the poloidal field line in the 
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(rO) plane. In a relativistic jet electric and magnetic fields are nearly equal in 
magnitude (t J2.1|) . E w | |, therefore the poloidal component of the energy 
flux is given by E\B^\I^k = (VLR/cfBl/Air. Applying this result at the BH 
horizon, the total electromagnetic power flowing out of a spinning BH is given 
by -fjet = k$1 ot £l 2 H , where 3>tot is the total magnetic flux threading the horizon, 
Qh is the angular frequency of the space time at the horizon, and k i s a constant 
which depends on the field geometry (see iTchekhovskov et alj|2009bl for details). 
Although this result for Pj e t is derived for a force-free jet in the limit of a 
slowly spinning BH (jBlandford fc Znaieklll977l ). it turns out to be surprisingly 
accu rate even for rapidly spin ning B Hs. and also for M HD jets (Tchekhovskoy 
et al. l2009bf ). 

One interesting question is whether the VL 2 H scaling of power is sufficient 
to explain the radio loud/quiet dichotomy of AGN. Observations indicate a 
factor of 1000 diff erence in radio power between radio loud and quiet AGN 
( Sikora et al.ll2007l ). Can this be entirely due to d i fference s in the BH spins of 
the underlying populations? ITchekhovskov et al. ( 2009b ) show that a model 
in which the BH is surrounded by a thick accretion disk, such that the jet 
subtends only a narrow solid angle around the axis, has a steep variation of 
jet power wi th spin - P < > 



et oc Vl H or even U° H (a similar steep dependence was 



observed by iMcKinnevl 2005 in the simulations of BHs with thick turbulent 
tori). Such a model is perfectly compatible with the observed radio loud/quiet 
dichotomy, whereas a BH with a thin disk around it would have P ]e t oc VL 2 H and 
is not likely to have a large range of jet power. This explanation of the radio 
loud/quiet dichotomy requires AGN jets, especially in low-luminosity systems, to 
be associated with thick advection-dominated accretion flows (ADAFs). There 



is, in fact, independent evidence for such a connection (jNaravan Sz McClintock 
20081 ). 



2.4 Stability 

According to the Kruskal-Shafronov criterion (e.g., Batemanlll978l ). cylindrical 



MHD configurations in which the toroidal field dominates are violently unstable 
to the m = 1 kink (or screw) instability. Since all models of relativistic mag- 
netized jets have B^ ^> B p , jets ought to b e highly unstable. This has been 



recently confirmed by numerical simulations ( Mizuno et al. 20091 ). Why are the 
jets observed in Nature coherent over eno rmous length scales? 

Narayan, Li, Sz Tchekhovskoy ( 20091 ) studied the stability of force- free jets 
and found that the growth rate of the kink mode is quite low. There are two 
reasons for this. First, although jets may have B^ 3> B p in the "lab frame," 
they typically have Ba ~ B p in the comoving frame of the fluid. This makes the 
instability less severe. Second, as a result of time dilation, the growth rate as 
measured in the lab frame is further suppressed by a factor of 1/7. Note that 
the above study was limited to a very simple force-free model. The stability of 
more realistic relativistic MHD jets is an open question. 



3 Relativistic MHD Jets 



We now consider MHD jets with gas inertia, though, for simplicity, we ignore 
gas pressure. A steady axisymmetric MHD flow has a number of conserved 
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Figure 1. [Panel (a)] Solid lines show Lorentz factor 7, magnetization pa- 
rameter a and conserved energy flux p, (= 450) as a function of distance 
r (roughly in units of the BH horizon radius) along an equatorial field line 
(8 = 90°) for a mass-loaded split-monopole outflow. Initially, 7 increases 
rapidly with r. However, beyond the fast point (indicated by the vertical dot- 
ted line), it grows only logarithmically. As a result, even at astrophysically 
relevant distances, 7 is only ~ 40 and most of the energy remains locked up 
in the magnetic field: tr ~ 10. This is the o problem. [Panel (b)] Similar to 
panel (a), but for a polar field line at 6 — 10°. The two vertical dotted lines 
show the positions of the fast point (left) and the causality point (right). Here 
the flow continues to accelerate beyond the fast point and up to the causality 
point. As a result, it is less affected by t he a problem and achieves a larger 
final Lorentz factor 7fi na i ~ p. (Based on iTch ckhovs kov et alj|2009af l 



quantities: (i) The angular velocity f2 is constant along any field line, (ii) The 
net mass flowing along a bundle of field lines is conserved, which means that the 
quantity "fpVp/Bp is constant along each field line, where p is the mass density 
in the rest frame of the gas and v p is the poloidal 3-velocity of the gas in the 
lab frame, (iii) The angular momentum flowing along the bundle of lines is 
conserved, (iv) Finally, the energy flowing along the bundle is also conserved, 
which implies that 

p = j|£? x B^/Airl + (7 2 /0V p ) j MpVp = constant, (1) 

where the first term in the numerator is the Poynting energy flux and the second 
is the gas kinetic energy flux. This equation may be rewritten as 

p = 7 (<T + 1), a=\ExB (f> /47T\/( 7 2 pv p ), (2) 

where the magnetization parameter a represents the ratio of electromagnetic to 
kinetic energy flux at any point in the flow. 

3.1 Acceleration and the a Problem 

Equation ([2]) has a simple interpretation. At the base of the jet, before the gas 
has accelerated, we have 7 ~ 1, and so the energy flow is almost entirely in 
the form of Poynting flux: a 3> 1, p S> 1. Here the jet is magnetically very 



Relativistic Jets 



5 



dominated and it behaves for all intents like a force-free jet. As the jet moves 
out, 7 increases and the kinetic energy flux grows at the expense of the Poynting 
flux. Thus o decreases as 7 increases, but in such a way as to keep fi constant. 
If the acceleration is efficient, then we expect a to become vanishingly small at 
large distance and all the initial energy to end up in the gas; we would then have 
7 = 7max = In fact, this ideal is rarely achieved. 

The problem is most clearly seen in the split mono pole geometry , where 
analytical results are available for equatorial field lines (Michel 1969; Beskin 



et al. Il998l ; iTchekhovskov et al.ll2009al ). The flow accelerates initially just as in 
the force- free case, with 7 ~ VtR/c. However, once the flow speed exceeds the 
fast magnetosonic wave speed in the medium, which happens when 7 ~ /U 1 / 3 , or 
equivalently when 7 ~ a 1 / 2 (from eq. [5] with /i S> 1), the gas loses contact with 
gas behind it and acceleration slows down drastically. Beyond this point, there 
is only a growth in 7 of order a logarithmic factor, so the final 7 is given by 
7final ~ Co 1 / 2 , where C is a logarithm. If /x is large, as it must be for a highly 
relativistic flow, the asymptotic value of a remains large and most of the energy 
is still carried in the form of Poynting flux rather than as gas kinetic energy. 
This is the a problem. 

The a problem is illustrated by a numerical example in Fig. [T^,. As we 
see, the problem is not a lack of energy. In this example, [x = 450 which means 
there is plenty of energy available and the gas could in principle accelerate up 
to 7 m ax ~ 450. However, it only accelerates up to 7fi na i — 40 because the energy 
remains stuck in the magnetic field. 

Simulations such as those shown in Figs. [TH2]of this paper are challenging 
and have only recently become possible with the use of the general relativis- 
tic MHD code HARM (Gammie et al. 2003), including recent improvements 



(jTchekhovskov et al. 2007, 2008). There are two difficulties with these calcu 



lations. First, the large value of a means that the electromagnetic and kinetic 
terms in the energy equation are very different in magnitude. The equations 
are thus stiff, making it very difficult to maintain accuracy. Second, magnetic 
acceleration is relatively slow: 7 varies only as z 1 / 2 even when acceleration is 
efficient, and it grows logarithmically when it is inefficient. Therefore, in order 
to obtain useful results, it is necessary to simulate relativistic jets over very large 
length scales, which is obviously challenging. As of this writing, only one other 
group (Komissarov and collaborators) is able to carry out such calculations. 

3.2 No a Problem for Collimated Jets! 

Observat ions sugges t that relativistic jets in Nature do not suffer from the a 
problem. Lind et al. ( 19891 ) showed that the interaction of a jet with an external 



medium is very different depending on whether <j < 1 or <r > 1. In the former 
case there is a clear hot spot (or working surface) where the jet shocks with the 
medium, and there is a well-defined backflow in a cocoon. This is similar to 
what is observed in FR II jets. In contrast, a magnetically dominated jet drills 
through the medium and has hardly any cocoon. It appears that AGN jets 
manage to achieve small values of a before they shock on the external medium. 
Why are they not limited by the a problem? 



A key clue was provided by ITchekhovskov et al.l ( 2009al ) who showed that 



equatorial and polar field lines in the split monopole geometry behave very dif- 
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Figure 2. The dotted lines in panel (a) show the variation of 7, u and 
fi (w 900) as a function of distance r (BH units) for an MHD jet that is 
confined by a rigid wall with a shape given by R oc z 5 / 8 . As for all externally 
collimated jets, we see that this model does not suffer from the a problem. 
Asymptotically, the jet achieves a large Lorentz factor 7 w /x and a low 
magnetization a <C 1. The dotted line in panel (b) shows the variation of 
the jet opening angle 9, and in panel (c) the product 76*. Note that 76 1 < 1, 
which is a general feature of fully collimated jets. Solid lines in the three 
panels correspond to a jet that is confined out to a distance r = 10 4 5 and is 
then allowed to move freely in vacuum. The jet undergoes a period of rapid 
acceleration just as it exits into vacuum (see panel a), but its opening angle 
does not change (panel b). As a result, asymptotically, the jet has nearly the 
same Lorentz factor as in the previous example, 7 ~ /i/2, a ~ 1, but now 
7# ~ 15 ( panel c). This model may b e relevant for understanding GRB jets. 
(Based on iTchekhovskov et alj[2009d ) 

ferently. This is illustrated in Fig. [TJa which shows that acceleration along a 
polar field line continues well past the fast magnetosonic point. In fact, acceler- 
ation slows down only after the gas crosses another critical point, the "causality 
point," beyond which it can no longer communicate with the axis. The intro- 
duction of this new critical point means that the final Lorentz factor of the gas 
is a function of the angle 9 between the poloidal field and the axis: 

7final sin0~C<7 1 / 2 , (3) 

where C is the same logarithmic factor as before. The smaller the value of 6, 
the larger the 7fi na i that the gas can achieve for a given value of fj,. Thus, flows 
that remain close to the axis accelerate more efficiently. 

This result is confirmed by numerical simulations of paraboloidal and other 
types of confined jets. The dotted lines in Fig. [2] show an example of a confined 
jet with fi ~ 10 3 . This jet accelerates with no trouble up to 7fi na i of nearly 10 3 , 
achieving an asymptotic a < 0.1. The lesson from this example is that, so long 
as a jet is provided adequate collimation, it will accelerate smoothly without 
being limited by the a problem. 
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3.3 How are 7 and 9 related? 

One additional interesting result is seen in Fig. [2j while 7 increases, the jet 
angle 9 simultaneously decreases. In fact, it appears that jO ~ 1, as noted by 
Komissarov et al. ( 20091 ). Do relativistic jets always satisfy this constraint? This 



can be tested with GRB jets, w here observations provide strong constraints on 
both 7 and 9 (e.g., iPiranlboOSh . 

Long duration GRBs typically accelerate to 7 ~ few hundred before decel- 
erating in their afterglow phase. At the same time, jet breaks observed in GRB 
afterglows indicate typical jet opening angles 9 ~ 0.05 radia n. Thus, GRB jets 
appar ently have 7$ ~ few tens. This appears to violate the iKomissarov et ah 



(2009) relation mentioned above, which is also confirmed by Fig. [2j 



Fortunately, there is a simple solution. The solid lines in Fig. [2] show 
a second simulation in which the jet is initially confined over a few decades in 
distance and it is then allowed to move freely in vacuum. This is meant to mimic 
a collapsar in which the jet is confined so long as it is inside the stellar envelope, 
but becomes free once it escapes from the star. In this numerical example, 
we see that the gas experiences a burst of acceleration just as it is deconfined, 
without changing its opening angle very much. Asymptotically, the jet satisfies 
the scaling given in equation Q with C ~ 10. The particular simulation shown 
here has 7fi na i ~ 500, #fi na i ~ 0.04, <7fj na i ~ 1, which is close to typical values 
observed in GRBs. (Note that <7fi na i is not directly measured. However, the 
fact that GRBs emit a substantial fraction of their enerrgy in prompt emission 
implies that their jets cannot be Poynting-dominated.) 

There is room for further work in this area. For instance, if observations by 
the Fermi Obser vatory routinely fin d 7fi na i > 10 3 ( Kumar Sz Barniol Duranll2009l ) 



and #fi na i > 0.1 ( Cenko et al.ll2009l ) for many GRBs, we would need to explain 



why these jets have 76? > 100. The constant C in equation (J3]) is a logarithmic 
factor and is unlikely to be much larger than about 10 — 20. Therefore, either we 
must accept that a S> 1, i.e., GRB jets are Poynting-dominated and somehow 
manage to convert a large fraction of their energy to promt gamma-rays, or that 
MHD is not the appropriate framework for understanding GRB jets. 



4 Conclusion 

The main message of the work described here is that collimation is the key 
to relativistic jets. If a jet is suitably confined by an external medium, it will 
accelerate without difficulty. While GRB jets have a natural collimating medium 
in the surrounding stellar envelope, the situation is less obvious for jets in AGN 
and XRBs. A disk wind is the most likely collimator in these systems, which 
suggests that many prominent jets are probably associated with geometrically 
thick accretion disks with strong winds and outflows, i.e., ADAFs of various 
kinds ( Naravan McClintock 20081 ). This connection is worth exploring further. 



Equation ([2]) shows that the maximum Lorentz factor that a jet can achieve 
is equal to /x. A highly relativistic jet requires a large value of fi, which means 
that the jet must start out very magnetically dominated at its base with negligi- 
ble mass-loading. There is very little understanding of how mass-loading works. 
Qualitatively, one imagines that field lines that are connected to the disk have 
a ready supply of disk gas for mass-loading, whereas lines that penetrate the 
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BH horizon are more likely to be mass-free. One thus suspects that the most 
extreme jets probably emerge from spinning BHs rather than from disks. 

Finally, in all of our discussion we assumed that a coherent magnetic field 
is already present in the system. The origin of this field is not understood. It 
could be advected in by the accretion disk from outside. However, the efficiency 
of such advection is poorly understood and the topic is controversial. 
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